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Context. Primordial black holes are unique probes of cosmology, general relativity, quantum gravity and non standard particle physics. 
They open a new window on the very small scales in the early Universe and also can be considered as the ultimate particle accelerator 
in their last (explosive) moments since they are supposed to reach, very briefly, the Planck temperature. 

Aims. Upper limits on the primordial black hole number density of mass = 5 10 14 g, the Hawking mass (born in the big-bang 
terminating their life presently), is determined comparing their predicted cumulative y-ray emission, galaxy-wise, to the one observed 
by the EGRET satellite, once corrected for non thermal y-ray background emission induced by cosmic ray protons and electrons 
interacting with light and matter in the Milky Way. 

Methods. A model with free gas emissivities is used to map the Galaxy in the 100 MeV photon range, where the peak of the primordial 
black hole emission is expected. The best gas emissivities and additional model parameters are obtained by fitting the EGRET data 
and are used to derive the maximum emission of the primordial black hole of the Hawking mass, assuming that they are distributed 
■ like the dark matter in the Galactic halo. 

Results. The bounds we obtain, depending on the dark matter distribution, extrapolated to the whole Universe (n /Jgtf (M Jr ) = 2.4 10~ 10 
to 2.6 10~ 9 are more stringent than the previous ones derived from extragalactic y-ray background and antiprotons fluxes, though less 
model dependent and based on more robust data. 

Conclusions. These new limits have interesting consequences on the theory of the formation of small structures in the Universe, since 
they are the only constraint on very small scale density fluctuations left by inflation. Significant improvements by data gathered by 
the FERMI y-ray satellite are expected in the near future. The interest of a generalisation of this work beyond the standard particle 
oo : model and in extradimensions is briefly alluded, 
-xf ■ 
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fl") 1. Introduction hereafter PBHs) is regarded as an inescapable consequence of 
Q\ ' cosmological models, specially inflationary onces (for a review 
q i Black holes are classically defined by the property that things, see Carr 20 05a). PBHs can thus serve as astrophysical probes of 
•• . including light, can go in but can't go out. However, this mere unification and their discovery, even if only minuscule ones are 
definition should be modified if quantum effects are taken into expected to show up through their radiation, would be a land- 
account, since they allow particle emission, with potentially ob- mark on the way t0 quantum gra vity. 
$Jj ■ servable effects (Hawking 1 1 9751 Page & Hawking |19761 l. The 
C3 ; Hawking radiation is in fact compatible with the well known Formed by pure gravity, PBHs require no special extension 
' result of general relativity that nothing can escape from inside of the standard model of particle physics and are predicted on 
the black hole horizon. Involving quantum fluctuations of the quite general grounds to be created in the very early Universe, 
vacuum it can be conceived as the creation of virtual particle- though the mechanism of their formation remains conjectural, 
antiparticle pairs that are disrupted by tidal forces in the vicin- If they exist, they are expected to be much more emissive than 
ity of the horizon. The particle (antiparticle) with positive total the massive ones. Those tiny black holes have not yet bee n dis- 
energy escapes to infinity while the antiparticle (particle) with covered in spite of considerable efforts (see Barrau et al. 120031 
negative total energy falls into the black hole. This hypothetical and Carr|2005b|for reviews), but they are still searched for since 
process points to profound connections between thermodynam- their discovery would confirm the Hawking theoretical construc- 
ics, quantum mechanics and general relativity. Unfortunately, as- tion. Even an upper limit on their number density is of inter- 
trophysical black holes, either stellar or supermassive, radiate so est since, being the first gravitational collapsed objects to form, 
weakly that they remain practically black. they could giv e precious information on the early Universe (Carr 
On the other hand, beyond the stellar black holes and su- 
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2005a, 2005b). 



permassive ones sitting at the centre of galaxies (Reid |2008), Particle emission from black holes has been predicted in a 

well motivated observationally, a very much fighter category of 4D spa cetime in the fram ework of semi-classical quantu m grav - 

pnmordial origin has been postulated (Carr |1974>. The forma- ity (Hawking [T975l Page [19761 [19771 Carr & McGibbon 1998 ). 

tion of black holes in the early Universe (primordial black holes, It has been rece ntly ex tended to extradimen sions in the brane 

scenario (Argyres et al. 1998 Emparan et al. 120001 Emparan & 

Send offprint requests to: R. Lehoucq Reall 2008] Kanti[2p04). We here focus on the 4D case, postpon- 
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ing the (4 + n)D case to a forthcoming paper. We only mention an 
interesting coincidence between 4D and 10 or 1 ID black holes. 

A first constraint has been imposed on the mean mass den- 
sity of M — M* primordial black holes (in units of the critical 
density, Q.pbh{M*)) through the extragalactic y-ray background 
integrated both over the PBH mass spectrum and the lifetime of 
the Universe (Pag e & Ha wking [19761 McGibbon & Carr [T99Tl 
Carr & McGibbon [19981 improved by Barrau et al. 120051 . This 
limit, however, depends on poorly known quantities. First, the 
initial mass function of primordial black holes is purely conjec- 
tural because the conditions in the early Universe before decou- 
pling are poorly constrained and large density contrasts could 
have arisen (Blais et al. 2003). Second, the true extragalactic 
gamma-ray background is difficult to extract both observation- 
nally, because it could be contaminated by very local inverse- 
Compton emission and instrumental backgrounds, and theoret- 
ically, due to unresolved AGN contribution. At last, a similar- 
ity between the intensities of the PBH galactic and extragalactic 
gamma-ray fluxes in this energy range has been advocated by 
Wright ( 1996 ), based on its analysis of the EGRET observations, 
which we confirm in this work. It can be readily be explained 
noting that the dark matter column densities toward the galactic 
centre (« 0.3 g/cm 3 ) is, by chance, of the same ordrer of mag- 
nitude than the extragalactic one integrated over a large part of 
the observable universe (« 0.2 g/cm 3 ). Thus it is reasonable to 
search for clues of primordial black holes in the diffuse galactic 
flux rather than in the extragalactic one, since it does not involve 
the redshift. 

A second constraint has been also obtained through the ob- 
servation of antiprotons in the galactic cosmic rays (Maurin et 
al. 120011 Barrau et al. 120021 and references therein). However, 
the limitation of his method is that the antiproton propagation is 
subject to uncertainties in their diffusion coefficient and the scale 
height of their confinement region. 

Here we search for the signature of PBHs in the Galactic y- 
ray background, more intense and better known than the extra- 
galactic one, taking benefit of the expected large concentration 
of PBHs in the central region of the Galaxy. More precisely, we 
assume that PBHs are broadly distributed like dark matter in the 
halo of our Galaxy. Admittedly, more sophisticated clustering as 
proposed by Chisholm (2006) should be analyzed in the future. 
This method requires a careful modelling of the Galactic back- 
ground, but it is free of the ambiguities of the initial PBH mass 
spectrum since it implies only black holes radiating now (i.e. of 
mass M* = 5 10 14 g for a 4D universe). It is not plagued either 
by the large uncertainty in the unresolved source contribution to 
the extragalactic background. 

2. Hypothesis and ingredients 

Our aim is to calculate the Galactic y-ray emission due to pu- 
tative PBHs assumed to be distributed as dark matter. We thus 
focus on those born during the very early phases of the Universe, 
13.7 Gyr ago, and dying presently. Their initial mass and temper- 
ature are estimated to M* = 5 10 14 g and T+ = 20 MeV (Carr & 
McGibbon [T9"9"8l MacGibbon & Webber [19^01 ). Because a black 
hole will emit all particle species with rest mass less than or 
equivalent to the black hole temperature, the emission is insensi- 
tive to the quark and gluon production which occurs above 200- 
300 MeV. The only species of interest in our case are light par- 
ticles (neutrinos, gravitons, photons and electron-positron pairs) 
if low mass non standard particles are excluded. Their relative 
proportions (branching ratios) are rather well established. 
The main hypotheses of this work are the following : 



1. spacetime is 4D; 

2. PBHs form through a cosmological scenario; 

3. most PBHs are presently neutral and non rotating; 

4. being part of the dark matter, PBHs are distributed alike. 

Like previous authors, we adopt points 1, 2 and 3 with- 
out discussion, point 4 is reasonable and generally assumed 
(MacGibbon & Carr [19911 Barrau 12002b . Two ingredients are 
crucial to perform the calculation: the y-ray emission spectrum 
of an individual black hole with mass M* and the Galactic dis- 
tribution of PBHs. The first one, relying on detailed calculations 
of McGibbon & Webber ( 119901 ) is rather well understood. The 
second one lies on a more shaky ground, since the distribution 
of PBHs in the Galaxy inherits all the uncertainties of the dark 
matter distribution. We will show, however, that this uncertainty 
does not invalidate our analysis. 

2. 1 . PBH emission spectrum 

The spectrum of the Hawking radiation emanating from a black 
hole with temperature 7h is not purely thermal since it has not an 
exact blackbody profile. The correcting factor, also called grey- 
body factor, depends on the energy of the emitted particles and 
strongly on their spin (high spin being defavoured). It signifi- 
cantly modifies the emitted spectrum below kTu. This distortion 
can be attributed to the fact that any particle emitted near the 
black hole horizon has to traverse a strong gravitational back- 
ground before reaching the observer at infinity, at variance with 
what happens with a black body in a flat spacetime (Page |1976l ). 
Consequently, the flux peaks at higher energy (around 5 kTu) 
than for a pure blackbody at same temperature (which flux is 
maximum at 1.59 kT); see FigureQ] The spin dependence of the 
greybody factor determines the relative emissivities (branching 
ratios) in different particle types from the black hole. 

Since the typical temperature of PBHs born in the early 
Universe and that end its life at present time is about 20 MeV, a 
distinctive signature of quantum black holes would be a quasi- 
planckian spectrum at unexpectedly high energy, peaking at 
about 100 MeV. This makes them particularly suitable to obser- 
vation by y-ray satellites. We thus reassess the question of lim- 
iting their maximal number in the light of EGRET observations 
of the diffuse y-ray emission of the Galaxy complementary to 
previous attempts relying on the extragalactic y-ray background 
(Page & Hawki ng [19761 Halzen et al. [19911 Carr & MacGibbon 
1998, Kim et al. [19991 Barra u et al . 120051 ) and antimatter in cos- 
mic rays (Boudoul & Barrau I2003I ). 

2.2. Dark matter density distribution 

It is reasonable to assume that PBHs are a minor part of dark 
matter and that, as such, they follow its calculated distribu- 
tion. Thus, the expected y-ray flux is proportional to the dark 
matter density which is quite different from the scenario based 
on neutralino annihilation where the y-ray flux varies quadrat- 
ically with dark matter density. The dark matter distribution, 
deduced from numerical simulations (Navarro, Frenk & White 
119971 Moore et al. 1999) is conjectural and quite uncertain. We 
use the most advertised dark matter distributions in the litera- 
ture. For the sake of simplicity, we do not discuss clumping in 
this article and consider only smoothed distributions. The inner 
Galaxy (« 100 pc from the center), where the dark matter density 
(and hence presumably the PBH density) is larger, is the main 
place for dark matter indirect searches through y-ray signatures. 
Indeed, since the total mass in the inner Galaxy is dominated 
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Fig. 1. Instantaneous y-ray photon emission rate from a kT = 
20 MeV black hole (McGibbon & Webber [T990t compared to 
the black-body emission at the same temperature (dashed line). 



by baryons, the dark matter distribution is likely to have been 
influenced by the baryonic potential (Bertone & Merrit 2005). 
Adiabatic compression models including the effect of the bary- 
onic gas significantly increase the dark matter density in the cen- 
tral region of the Milky Way, and as such, lead to stronger con- 
straints on the PBH number density. More realistic than the stan- 
dard model, they are unfortunately still quite uncertain. The halo 
models with adiabatic compression have been discussed and put 
under a useful form by Mambrini et al. ( 120061 see their table 
1), from which we borrow the results. The Navarro-Frenk- White 
(hereafter NFW) distribution diverges in the Galactic centre di- 
rection. To cure this divergence we use the prescription from 
Delahaye et al. d2008l 

Note that a 4 millions solar masses supermassive black hole 
sits at the centre of the Milky Way. Adiabatic accretion of dark 
matter onto this object could lead to a very high density of dark 
matter (Bertone & Merrit 2005 ). An initial dark matter density 
profile p{r) oc r~ a would lead to a final density profile with a 
slope a' — (9 - 2a/ (4 - a) so sharp that it is equivalent to a point 
source for EGRET and Fermi satellites. 



3. Modelling the y-ray emission from the interstellar 
medium and PBHs in the Milky Way 

Our aim is to estimate the maximum y-ray flux from PBH emis- 
sion in the Galaxy that is allowed by observations on top of 
the cosmic-ray induced y-rays of interstellar origin. More pre- 
cisely, we compare the total intensity expected from the vari- 
ous PBH distributions and known interstellar components to the 
whole-sky survey data obtained by the Energetic Gamma-Ray 
Experiment Telescope (EGRET). The latter, which operated on- 
board the Compton-Gamma Ray Observatory from April 1991 
to May 2000, detected photons in the 20 MeV to 30 GeV range. 
The observation program made use of the large instrumental 
field of view to cover the entire sky and to complete in-depth 
studies of specific regions. EGRET is well suited for this study 
since the peak energy of PBH emission (E » 100 MeV) is inside 
the range of high effective area of the detector (Thompson l 1 9931 1. 
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Fig. 2. Dark matter density profiles following NFW dl997l l, 
Moore et al. dl999b without and with adiabatic compression 
(NFW C and Moore c respectively, Mambrini et al. 2006). An 
isothermal profile is also added for comparison. All these distri- 
butions are normalized to 0.3 GeV/cm 3 in the solar neighbour- 
hood. 



Moreover, EGRET repeatedly observed the Galactic centre and 
Galactic plane regions which thus benefit from long exposures 
and limited off-axis aberrations. 

The Galactic emission is produced by the interaction of en- 
ergetic cosmic-ray electrons and protons with interstellar nucle- 
ons and photons. The decay of neutral pions produced in hadron 
collisions, the inverse-Compton scattering of the interstellar ra- 
diation field by electrons and their bremsstrahlung emission in 
the interstellar gas are the main contributors to the Galactic 
emission. If energetic cosmic rays penetrate uniformly all gas 
phases, the y-ray intensity in each direction of the sky can be 
modelled by a linear combination of gas column-densities, an 
inverse-Compton intensity map (lie), an d an isotropic intensity 
(liso) which accounts for the very local inverse-Compton emis- 
sion, the extragalactic background, and the residual instrumen- 
tal background. To account for the non-uniform cosmic -ray flux 
in the Galaxy, the gas column densities are distributed within 
galactocentric rings. The all-sky Leiden- Argentina-Bonn (LAB) 
composite survey (Kalberla et al. 2005) was used to obtain the 
atomic hydrogen Nhi column densities assuming a constant spin 
temperature of 125 K. The molecular hydrogen column den- 
sities were inferred from the velocity integrated CO intensity 
W(CO) obtained from the Center for Astrophysics compilation 
(Dame et al. |2001| l. Following Strong and Mattox ( 119961 ), we as- 
sumed a uniform Xqo = N{Hz)/W(CO) conversion ratio with 
galactocentric distance to reduce the number of free parame- 
ters in the model. The gas column densities within 10° of the 
Galactic centre and anticentre were interpolated from the re- 
gions just outside these boundaries, then normalized to match 
the total emission observed along the line-of-sight. We have also 
included within the solar-circle ring the large column-densities 
Nudark of "dark" gas associated with cold and anomalous dust 
at the transition between the atomic and molecular phases in the 
nearby clouds (Grenier et al. 2005). To derive the spatial distri- 
bution and intensity of the inverse-Compton radiation, we have 
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used the GALPROP model for cosmic -ray propagation devel- 
oped by Strong et al. ( 120071 1, more precisely the version 50p 
of this model with an electron source spectrum similar to that 
measured near Earth. This version assumes a cosmic-ray source 
distribution matching the radial distribution of Galactic pulsars 
from Johnston (1994). Maps of PBH emission, Ipgfiihb), were 
constructed as a function of energy using the spectrum displayed 
in FigureQ]for the various radial profiles of dark matter densities 
given in Figure [2] Our model also includes the flux and spatial 
distribution of all point sources that were detected above 6<x in 
the revised EGRET catalogue (Casandjian & Grenier 2008 ) . 

The predicted count rates per bin in a model map are there- 
fore calculated as: 

NpredQ, b) = 6(7, b) [qichciU b) + /„„ + qdarkNudarkd, b) 

+ Zi=rm g s qmANmin, I, b) + 2X C0 W co (n, I, b)] 
+qpBHlpBH(l, by} + 2 j-sources 

e(lj,bj)fjPSF(lj,bj) 

Here fj and e(l, b) represent source fluxes and the EGRET 
exposure as a function of Galactic coordinates, respectively. The 
product of the diffuse maps with exposure were convolved with 
the EGRET point-spread function (PSF) for an input E~ 2 1 spec- 
trum before adding the source maps. The EGRET data, expo- 
sure, and instrument response functions were downloaded from 
the CGRO Science Support Center. The q parameters (gas emis- 
sivities or relative contributions of different radiation compo- 
nents), the isotropic intensity, and the Xqo factor were fitted 
to the EGRET photon maps in six energy bands, namely 70- 
100 MeV, 100-150 MeV, 150-300 MeV, 300-500 MeV, 500- 
1000 MeV, and 1000-2000 MeV. We used a 2-dimensional 
maximum-likelihood method with Poisson statistics using maps 
in 0.5° x 0.5° bins. The likelihood (L) is calculated as the prod- 
uct, for all pixels at 10° < |/| < 170° and all latitudes, of the 
Poisson probabilities of observing N bAU b) photons in a pixel 
where the model predicts N pre d(l, b). 

Figure [3] shows the EGRET count map in the merged 70 
to 150 MeV band, together with the predicted count map ob- 
tained for the same energy range from the best-fit model with- 
out any PBH contribution. The bottom map gives the residuals 
expressed in cr units, (N b s (l, b) - N pre d{k b))/N pre d(l, b) 1/2 . The 
latter is smooth and shows no obvious positive residual toward 
the inner Galaxy, so the combination of interstellar, isotropic, 
and point-source components reasonably fit the observed y-rays. 
The patches of small, but systematic, residuals spanning angular 
scales of 20 to 30 degrees in various places across the sky are 
correlated with the exposure map and are due to off-axis aber- 
rations and systematic uncertainties in the off-axis effective area 
of the telescope. 

A useful constraint on the PBH number density can be de- 
rived if their cumulative emission in the Milky Way is well sep- 
arated in space and energy from the main Galactic and isotropic 
contributions. The prominent lens-like source of y-rays seen in 
the data and model maps at the centre of Figure [3] results from 
the enhanced cosmic -ray density and interstellar radiation field 
in the inner Galaxy. The photon field includes the stellar light 
from the Galactic bulge which is concentrated within about 10 
degrees of the Galactic centre (Moskalenko et al. 2006 ). Because 
dark matter piles up near the centre as well, one expects some 
degeneracy between the integrated intensity from PBHs and the 
inverse-Compton radiation, but they exhibit different ellipticities 
(latitude to longitude ratios) that allow their separation in the fit 
(see Figure @}. The intensities in Figure [4] are integrated from 
70 to 150 MeV where the PBH spectrum peaks, but where the 
EGRET angular resolution is poorer. It illustrates that the two 
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Fig. 3. From top to bottom: EGRET photon map, best-fit pho- 
ton map and residuals. The first two maps are expressed in pho- 
tons/bin while the latter is in cr units. All maps are given in 
Galactic coordinates for the 70 to 150 MeV energy band. The 
data and model maps are displayed in logarithmic scale, the 
residuals in linear scale. 



maps still differ after convolution with the wide instrument func- 
tion. Figure [5] shows the relative longitude profiles, after con- 
volution with the PSF, of the inverse-Compton component and 
PBH emission for various dark matter distributions. All PBH 
profiles but the isothermal one have smaller angular scales than 
the inverse-Compton one. The rapid drop of the isothermal dis- 
tribution outside 30 degrees, together with a different spectrum, 
reduces the degeneracy between the two components. The PSF 
profile, which is also displayed in Figure [4] appears to be sig- 
nificantly sharper than even the most concentrated dark matter 
distributions, thus allowing separation of the PBH contribution 
from a soft point-source. 

The likelihood probability can serve to test the presence or 
absence of PBH emission in addition to the other radiative com- 
ponents. We derived the likelihood ratio test statistic defined as 
TS = 2[ln(Li) - ln(Lo)L where the likelihood values L\ and Lq 
are optimized respectively with and without a PBH contribution 
in equation [3] Asymptotically, the TS distribution follows &x\ 
distribution and the detection significance of PBH emission is 

VTS" cr. 

The y/TS results obtained in the six energy bands are given 
in Figure [6] They yield no significant detection at the 3 cr level 
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Fig. 4. Predicted photon counts in the merged 70 to 150 MeV 
energy band originating from 7r°-decay and bremsstrahlung radi- 
ation in the HI and H2 gas (top), inverse Compton scattering of 
the interstellar radiation field (middle), and PBH emission from 
the Moore dark matter distribution (bottom). All maps are dis- 
played in Galactic coordinates with different logarithmic scales 
to easily compare the radial profiles. 
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Fig. 5. Normalized longitude profile of the EGRET PSF (black 
continuous thick), exposure (black dotted) and inverse-Compton 
(IC, black dashed) maps, as well as the profiles of the PBH 
emission maps obtained for the dark matter distributions of 
Moore (red), Moore c (green), NFW (black thin), NFW C (blue), 
and isothermal (magenta). All profiles are given for energies 
E > 100 MeV and are centered on b = 0°. The IC and PBH 
profiles include the convolution by the PSF at these energies. 
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Fig. 6. Comparison of Galactic emission model, with and with- 
out the PBH contribution, to the EGRET count map in six energy 
bands. The result expressed in yTS, or detection significance in 
unit of cr, is shown for five dark matter distributions: Moore, 
Moore c , NFW, NFW C and isothermal (from top to bottom). 



of any PBH emission for the dark matter distributions of Moore, 
Moore c , NFW, and NFW C distributions. A small likelihood in- 
crease, amounting to 6 cr, is found for the isothermal distribu- 
tion at low energy. It vanishes above 300 MeV. Although a low- 
energy excess is expected from black hole emission, the signal 
is not significant enough to be interpreted as a signature of PBH 
activity. 

The main reason is that the intensity does not decrease in 
the 150-300 MeV energy range as it should (it should amount 
to only 7.5% of the total PBH intensity). We also point out that 
the isothermal distribution has the nearest angular profile to the 
inverse-Compton one. The additional emission that is attributed 
by the fit to the isothermal map is likely to compensate for a lack 
of inverse-Compton emission in the central part of the Galaxy. 
One should then note that gas in the dense central molecular 
zone must be mapped using rarer CO isotopes which are not 
part of the present gas maps. The central molecular zone (CMZ, 
about 500 x 500 x 50 pc) has a high average hydrogen den- 
sity (n H * 50 to 100 cm" 3 , Launhardt et al l20UZl Femere 120071 
about 90% of the interstellar gas being in the form of dense 
clouds (10 3 cirT 3 ) with higher values in cloud cores. The remain- 
ing 10% would form a diffuse homogeneous intercloud medium 
easily pervaded by GeV cosmic rays (which produce most of 
the 100 MeV photons). The Galactic centre region also abounds 
in massive stars, in particular in young and dramatically active 
OB associations (the Arches, Quintuplet, and Central clusters, 
Figer et al. 120001 |2008l > where SNII shocks and colliding stellar 
winds from densely packed Wolf-Rayet and O stars are favorable 
to cosmic-ray acceleration. So, we can expect both a strong en- 
hancement in the y-ray intensity corresponding to CMZ and its 
surroundings when cosmic rays diffuse out, and a different spec- 
trum from freshly accelerated cosmic rays compared to the local 
one. Both these effects may have been observed at TeV energies 
and in the radio (see Grenier 2008] for a review), but they have 
not been included in the present model. Furthermore, the 2<x 
likelihood increase found above 1 GeV for all PBH distributions 
may be related to the present use of the Earth in-situ cosmic- 
ray electron spectrum to derive the inverse-Compton intensity, 
whereas the EGRET data requires a harder IC spectrum to fit the 
data. This GeV excess, probably of instrumental origin (Stecker 
et al. |2008l l, is not taken into account in the present model. Given 
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these significant uncertainties in the actual shape and spectrum 
of the inverse-Compton emission in the inner Galaxy and from 
the presence of an additional contribution from the CMZ gas, we 
do not consider the likelihood increase related to the isothermal 
PBH component as significant. 

To obtain a fair upper limit on PBH emission, we merged the 
70-100 MeV and 100-150 MeV energy bands where 72% of the 
signal is expected. We used the likelihood test ratio to increase 
the PBH contribution until the Log-likelihood drops by 9 from 
the maximum value obtained under the null hypothesis (without 
PBH emission). As required by the definition of TS, all other free 
parameters (interstellar, isotropic, and sources) were optimized 
as we changed the PBH intensity. In other words, we increased 
the PBH contribution until the predicted map was so degraded 
that there was a 10 3 chance that random fluctuations from the 
null+PBH model yield as good a fit to the data as the null model. 
The derived upper limits are presented in Table [TJ 

4. Results 

We have determined the maximum fractions, /, of dark matter 
in the form of PBHs in our Galaxy which decay presently and 
participate to the Galactic y-ray emission for different possible 
dark matter distributions (see Table [TJ. Knowing the dark mat- 
ter density in the Universe, Q^m = 0.228 ± 0.013 (Komatsu 
et al. 120081 1. these fractions can be translated into Q.pbh{M*), 
the cosmological density of PBHs with mass equal to M+: 
Qpbh(M*) ~ f^DM- We get values ranging from 2.4 10~ ll) to 
2.65 10~ 9 which leads to a significant improvement for the con- 
centrated Moore c and NFW C distributions over previous bounds 
deduced from the extragalactic y-ray background, i.e. 10~ 8 (Carr 
& McGibbon [T998l and 3 10~ 9 (Barrau et al.|2003]) and f rom the 
number of antiprotons, 3 to 4 10~ 9 (Maurin et al. 120011 Barrau 
et al. 120021 ) . It is worth noting that our method is independent 
of the poorly known initial mass function of PBHs. Our results 
are consistent with a completely independent constraint derived 
from the apparent absence of gravitational capture of PBHs with 
M > M+ by solar-type stars and white dwarfs (Roncaldi et al. 
120091 . 



Table 1. Constraints on PBHs of mass M* for different dark mat- 
ter distributions 3 . 



DM distribution 








Moore 


6.04 ± 0.05 lO" 9 


1.38 10-" 


0.98 10- 2 ' 


Moore c 


1.07 ± 0.07 10~ 9 


0.24 10~ 9 


0.17 10~ 27 


NFW 


6.70 ± 0.05 10~ 9 


1.53 10~ 9 


1.08 10~ 27 


NFW C 


1.93 ± 0.08 10-' 


0.44 10- 9 


0.31 10- 27 


isothermal 


11.62 ±0.04 10-" 


2.65 10~ 9 


1.87 10~ 27 



a / is the maximum fraction of dark matter in the form of PBHs, 
£Ipbh(Mj,) is the corresponding cosmological density, and j} is the 
fraction of regions undergoing collapse when the mass enclosed 
within the cosmological horizon was M t . 



5. Astrophysical consequences 

5.1. Local PBH density 

We can derive directly the local density of PBH of mass (and 
the mean distance between them) from the dark matter density 
in the solar vicinity (i.e. pdm(O) ~ 0.3 GeV/cm 3 ). We get values 



ranging from 3.2 10~ 1() to 2.0 10~ 9 GeV cirT 3 corresponding to a 
maximum number density of M* PBHs within the range 3.3 10 7 
to 2.1 10 8 per pc 3 . These values, in turn, can be used to derive an 
upper limit of the present explosion rate of PBHs (Schroedter et 
al. 120081 Petkov l2U08l and references therein). 

To estimate the local explosion rate of PBHs, we rely on 
the Wright's pioneering analysis of the y-ray PBH emission of 
the dark halo (Wright 1996). This author deduces that the con- 
centration factor of PBH, which is the ratio of the local to the 
cosmic PBH densities should be less than 2 - YlhT x 10 s pc -3 , 
i.e.2.74 10 5 - 1.64 10 6 pc~ 3 taking the reduced Hubble con- 
stant h = 0.73. Combining this estimate with the work of 
Halzen et al. ( 119911 ), Wright deduces a limit on the explosion 
rate < 0.07 - 0.42 pc~ 3 yr _1 . Since our concentration factor, 
is =s 2.3 10 5 (because Pdm(Q) ~ 0.3 GeV/cm 3 while pdm ~ 
0.23 p c 1.3 10~ 6 GeV cm -3 , where p c is the critical density), 
by simple scaling we get an explosion rate 0.059 pc~ 3 yr _1 . 

This limit should be compared to the observational ones. 
However the behaviour of high temperature PBH (which has low 
mass) is poorly known, and different models of the final burst 
can give vastly different detection limit. The final fate of PBH 
and the ability to observe their final decay depends on the model 
of particle physics chosen to describe their high temperature be- 
haviour. In the standard model the number of elementary particle 
species never exceeds 100, and the likelihood of detecting such 
explosions is very low. However, the physics of the QCD phase 
transition around 157-180 MeV is still uncertain. The prospect 
of detecting PBH explosions would be improved in less con- 
ventional particle physics models. For instance in the Hagedorn 
model the number of particle species, including resonances, pro- 
liferates at the QCD temperature. The resulting burst of radiation 
is therefore much shorter and intense than in the standard case. 
In the standard model case, the limit derived here on the local 
rate density of exploding PBH provides a very difficult target 
for all techniques to direct detection of PBH explosions through 
y-rays. By comparison, our explosion rate is: 

- near the EGRET observational limit (0.05 pc~ 3 yr _1 , see 
Fichtel et al. 119941 ), assuming that the last 6 10 13 g of PBH 
rest mass evaporates, producing 10 34 ergs of 100 MeV y-rays 
in less than a microsecond, based on the Hagedorn model of 
high energy, which remains the most stringent y-ray obser- 
vational limit at present time; 

- higher than the Fermi expected detection rate 
(» 0.02 pc~ 3 yr _1 for 2 years observations, see Schroedter et 
al. 120081 : 

- higher than the Maki et al. 119961 ) antiproton limit 
(0.017 pc~ 3 yr~') but independent of the delicate cosmic ray 
diffusion coefficient, since y-rays travel in straight lines. 

5.2. Multidimensionals PBHs 

Interestingly enough, extending this study to higher dimensions 
(Argyres et al. 1998), we find that the temperature of the PBH in 
a 11D spacetime, with a lifetime equal to that of the Universe, 
is close to that of 4D ones. Indeed their Schwarzschild radius 
(5 fm) is very close to the radius of the extradimensions (7 fm) 
if the fundamental Planck mass is of order 1 TeV, as suggested 
by Arkani-Hamed, Dimopoulos and Dvali ( 1998). Consequently, 
they are expected to behave like "normal" 4D black holes. 
For such a high number of extra dimensions, the Kaluza-Klein 
modes of the graviton are so energetic (being proportional to the 
inverse radius of the extra dimensions) that they are not excited 
at the PBH temperature. Thus the branching ratios of a 4D black 
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hole can be applied to them safely. Their calculated temperature 
being 24 MeV they are almost undistinguishable from 4D PBH. 
More work is necessary to settle the question since when the 
black hole radius is of the order of the Planck length we enter in 
the domain of pure quantum gravity, poorly understood. Thus, 
qualitatively at least, the conclusions of this work could apply to 
superstring black holes. We will discuss this issue in much more 
details in a forthcoming paper. 

5.3. Cosmological consequences 

PBHs should be the first objects to form by gravitational col- 
lapse of density inhomogeneities in the early Universe, and as 
such they are a unique tool to probe the very small cosmolog- 
ical scales, an order of magnitude smaller than that explored 
by measurements of the cosmic microwave background and/or 
large scale s tructur es (Carr & MacGibbon ll998l Khlopov 2008, 
Barrau et al. 120031) . The fraction of regions undergoing collapse 
when the mass enclosed within the cosmological horizon is M 
is determined by the equation of state of the medium in which 
PBHs form (P = yp) and the root mean square amplitude e 
of the fluctuations entering the horizon at that epoch (Carr & 
MacGibbon l 19981 ) through 

y8(Af) = e(M) exp(-y 2 /2e(M) 2 ), 

provided the fluctuations are gaussian and spherically symmetric 
and PBH production takes place in the radiation dominated era 
(y = 1/3). Thus (i(M) offers a link between Q. PBH (M) and e(M). 
Numerically, 

CWM) * lO'^M) U^-J . 

Several constraints on (3 have been derived in different mass 
ranges from primordial nucleosynthesis and cosmic microwave 
background, i.e. for PBHs of mass smaller than (see Figure 1 
of Can- & MacGibbon |1998l l. Our values of Q/> B// (M*) translates 
to /3(M*) values listed in Table Q] Blais et al. ([2003 ) propose a 
different Q.-/3 relation based on a model in which the power spec- 
trum at small scale is boosted and where 

CW(M)/i 2 = 6.35 10 16 /3W (j^y ' ■ 

Our upper limit associated with h = 0.7 (Komatsu et al. 2008 ) 
leads to /3(M*) values around 10 times smaller than in the pre- 
vious case. One natural source of fluctuations would be inflation 
and, in this context, e(M) depends implicitly on the inflation po- 
tential, but this subject is beyond the scope of the present work. 

6. Conclusion and perspectives 

Our constraint on the maximum density of primordial black 
holes is more stringent, or at least comparable, than the previous 
ones. The method proposed is more direct and robust due to the 
fact that Galactic y-ray astronomy around 100 MeV has reached 
a level of unprecedented quality both from the theoretical and 
observational points of view. 

On the observational side, the large improvement in both 
sensitivity and angular resolution of the new generation Large 
Area Telescope (LAT) on board the Fermi Gamma-ray Space 
Telescope will allow to strengthen the upper limit on Q.pbh(M*). 
The central regions of the Galaxy will remain difficult to model 
precisely, but the sharper PSF will help separate low-energy 



PBH emission from the wider scale of the inverse-Compton ra- 
diation and from a population of point-sources in the Galactic 
bulge. The low-energy threshold of the LAT (20 MeV), the large 
photon statistics and very smooth exposure resulting from the 
continuous all-sky survey currently underway will help take ad- 
vantage of the unusual spectrum of PBH emission to separate it 
from the harder and power-law spectra of inverse-Compton and 
gas emission, or from the more isotropic contamination from soft 
instrumental backgrounds. The low statistics of EGRET has not 
permitted to fully exploit this aspect. 

On the theoretical side, the extension of this method to non 
standard light particles would shed new light on particle physics. 
Indeed, a PBH radiates any and all particles with rest masses 
which are substantially less than its current temperature, includ- 
ing graviton and hypothetical light dark matter particles (Boehm 
et al. 120041 Fayet 120041 ). axions and light scalars. Concerning 
light fermionic dark matter invoked to explain the low-energy 
positron excess deduced from the 511 keV line observation of 
the inner Galaxy (Boehm et al. 120041 Fayet 120041 ). it would 
correspond to introducing an additional neutrino species which, 
opening a new decay chain, would slightly lower the y-ray yield 
and alter the lifetime of PBHs. The effect of non-standard par- 
ticles on the black hole radiation is postponed to a forthcoming 
paper. 
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